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Abstract 



^ I In this article we exhibit exphcitly the matrix model {6 = c«) fixed point of phi- 

^n I four theory on noncommutative spacetime with only two noncommuting directions using 

^f^ ' the Wilson renormalization group recursion formula and the 1/A^ expansion of the zero 

C^^ ■ dimensional reduction and then calculate the mass critical exponent u and the anomalous 

o 
en 



■^ ■ dimension ri in various dimensions . 

o ' 



1 Introduction 



The Wilson recursion formula is the oldest, most simple and most intuitive renormaliza- 
^_- tion group approach which although approximate agrees very well with high temperature 

expansions [1]. 

The goal in this article is to apply this method to the self-dual noncommutative phi-four 
in the matrix basis [3] which after appropriate non perturbative definition becomes an A'^ x 
N matrix model where A^ is a regulator in the noncommutative directions. More precisely 
we propose to employ, following [2l[10l[TT], a combination of i) the Wilson approximate 
renormalization group recursion formula and ii) the solution to the corresponding zero 
dimensional large N counting problem given in our case by the hermitian Penner matrix 
model which can be turned into a multi trace hermitian quartic matrix model for large 
values oi 9. As discussed neatly in [2] the virtue and power of combining these two 
methods lies in the crucial fact that all leading Feynman diagrams in 1/A^ will be counted 
correctly in this scheme including the so-called "setting sun" diagrams. As it turns out 
the recursion formula can also be integrated explicitly in the large N limit which in itself 
is a very desirable property. 
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In a previous work |18j a non perturbative study of the Ising universality class fixed 
point in noncommutative 0{N) model was carried out using precisely a combination of 
the above two methods. It was found that the Wilson-Fisher fixed point makes good sense 
only for sufficiently small values of up to a certain maximal noncommutativity. In the 
current work we focus on the opposite limit of large 9 although the 1/A^ expansion invoked 
in this article is different from the 1/A^ expansion of the 0{N) vector model since N here 
has direct connection with noncommutativity itself. The central aim of this article as we 
will see is to exhibit as explicitly as possible the matrix model fixed point which describes 
the transition from the one-cut (disordered) phase to the two-cut (non-uniform ordered, 
stripe) phase in the same way that the Wilson-Fisher fixed point describes the transition 
from the disordered phase to the uniform ordered phase. 

This article is organized as follows. In section 2 we write down noncommutative phi- 
four in the matrix basis and discuss some useful approximations involving the propagator 
at the self-dual point which are valid for large 0. In section 3 we consider the dimensional 
reduction of noncommutative phi-four and some of its properties. In particular we will 
derive a very simple non perturbative equation for the 2— point proper vertex which will be 
used to test the results obtained later using the 1/N expansion and the recursion formula. 
In section 4 we perform the tedious task of deriving the renormalization group equations 
which control the flow of the coupling constants of the model (three in this case) using the 
Wilson renormalization group recursion formula and 1/A^ expansion. The most difficult 
piece of the calculation as we will see is wave function renormalization. In section 5 we 
derive the non trivial fixed point and the associated critical exponents u and rj by solving 
numerically via the Newton-Raphson algorithm the renormalization group equations and 
discuss some of relevant physics. We conclude in section 6 by a summary of the obtained 
results and brief outlook. 

2 The Noncommutative Phi-Four Theory 

Let us consider a phi-four theory on a noncommutative Moyal-Weyl space with only 
two noncommuting coordinates, viz Rg = R x R^ where D = d — 2. We also introduce 
non-commutativity in momentum space by introducing a minimal coupling to a constant 
background magnetic field Bij as was done originally by Langmann, Szabo and Zarembo in 
[3] . This corresponds to the addition of a harmonic oscillator potential to the kinetic action 
which modifies and thus allows us to control the IR behavior of the theory. A particular 
version of this theory was shown to be renormalizable by Grosse and Wulkenhaar in [1] . 
The action of interest is given by 



J [ \ 2*2^^-'^ 2^2/ 4! 



(2.1) 



The index i runs over the noncommuting directions while the index fi runs over the 
commuting directions. We rewrite this action as 



S = d'^x 



2 



(2.2) 



The harmonic oscillator coupling constant Q, is defined by O^ = B'^O'^/A whereas the 
coordinate Xj is defined by Xi = 2{6^^)ijXj. It was shown in [5] that this action is 
covariant under a duality transformation which exchanges among other things positions 
and momenta as pi -^ xi. The value O^ = 1 in particular gives an action which is invariant 
under this duality transformation. Under the field/operator Weyl map we can rewrite the 
action as 



S = u,jd^xTrn[^^^-\d^ + \^^X^-\dl + ^y + ^^^^^^^ 



(2.3) 



The Planck volume V2 is defined by V2 = \/det27r0. We can expand the scalar fields in the 
Landau basis {(t>m,n} as (with x standing for commuting coordinates) 

^ oo ^ oo 

^ = ^= 5^ M„„(x)0„,„ ,$+ = —= Y. ^*mn^t^,n- (2-4) 



^2 '-^ , V^2 , 

m,n=l * 771,71=1 



The Landau basis is constructed for example in ^. The infinite dimensional matrix 
M should be thought of as a compact operator acting on the separable Hilbert space 
H = 5(N) of Schwartz sequences {am)m>i with sufficiently rapid decrease as m — > oo. 
In the Landau basis the action becomes 



j£> 



d xTru 



-d^M+d^M + -fi'^M+M + -r'^E{M, M+} + —{M+Mf + remainder . 

(2.5) 



The coupling constants r^ and u are defined by 

:^2 + i) 

V2 4! V2 

The remainder is given by 



2 47r(0^ + 1) \N 

r = ,u = - — . (2.6) 



1 1 -O^ 
remainder = -r'^^/IJ(r+M+rM + M+T+MT) , ^ = -^ . (2.7) 

The matrices F and E are given by 

\ 

{V)lm = Vm - 16lm-l , {E)im = {I - -^)^lm- (2-8) 

We regularize the theory by taking Af to be an A^ x A^ matrix. We will also restrict the 
momenta along the commuting directions to be less or equal than a hard cutoff A. The 
states (pi^m{x) with 1,171 < N where N is some large integer correspond to a cut-off in 
position and momentum spaces [1]. The infrared cut-off is found to be proportional to 
R = V29N while the UV cut-off is found to be proportional to An = y^8N/9. In [3j a 
double scaling strong noncommutativity limit in which N/9 (and thus An) is kept fixed 
was considered. In our analysis here we will take a different strong noncommutativity 



limit in which we take 9 — > oo first and then A^ — > oo. The trace Tr^ will be replaced 
by Tr]\f with Trj^l = N . 

This is a deformed principal chiral model. Indeed if we set the remainder equal to 
zero and choose /i^ = —Au/N we see immediately that in the limit u — > oo the partition 
function is localized around configurations satisfying M~^M = 1 with an action given 
by {J d^x Triqdp^M^dp^M)/2. This is the usual definition of a principal chiral field in 
dimension D = d — 2. The potential term provides a smoother definition of the model 
while the remainder gives a deformation due to non-commutativity. 

The free propagator simplifies in the limit 6 — > oo. Indeed in this limit the above 
action reduces to a U{N) matrix model in D = d — 2 dimensions with propagator given 
by l/(p^ + /i^) [TJ. A more interesting limit is Q,^ — > 1. This is the self-dual theory in 
which the off diagonal quadratic terms vanish, i.e. the remainder is 0. The action in this 
case is given by 



^D 



d xTrj\f 



2 ' ^2 2 ' ' N 



d^,M+d^,M + -^?M+M + -r'^E{M,M+] + —{M+Mf 



{2. 



We will consider in the remainder only the case of hermitian matrices M = M^ . There 
are three independent parameters in this theory. These are the usual mass parameter /x^ 
and the quartic coupling constant u plus the inverse noncommutativity r^ = 4/0. The 
free propagator of this theory is also simple given by 

/d^v 
-——^A,nn{p)e'P^''~y\ (2.10) 



p^ + /U^ + r^(?n, + n — Ij 

In the next section we will need extensively the sum Ylm^"^n.(p)- For 9 = oo this sum 
is obviously of order A^. Including also the subleading 1/9 correction this sum takes then 
the form 

Y^ A^jip) -^ NAnojip) , no = ^^. (2.12) 

m 

A straightforward generalization of (I2.12|] is 

5^ A-^,(pi)A-.^(p2)... -^ iVA-,,(pi)A-^.^(p2).... (2.13) 

m 

The structure of the 1/A^ expansion is quite complicated for generic values of 9. However, 
in the limit in which we take 9 — > oo first and then N — > oo we find that the 1/A^ 
expansion becomes manageable. More precisely if we assume that 9 is sufficiently large 
so that we are allowed to use the approximations (I2.12p and (I2.13P then the Feynman 
diagrams that will dominate the large N limit are exactly those of the hermitian matrix 
model in D dimensions. 



3 Quartic and Penner Matrix Models 

The reduction of the model ()2.9p to zero dimension yields the pure matrix model (with 
VD the volume of momentum space) 



V = voTrN 



-f/M^M + -r'^E{M, M+} + '^{M^Mf 



(3.1) 



We will be interested in the real hermitian matrix model with M+ = M. By an appro- 
priate rescaling we rewrite this as 



Vm = Ttn 
The coupling constants m? and g are given by 



-M^ + m^EM'' + j-M'' 



(3.2) 



2 r^ vdu 



The path integral is given by 



"^ =— ' 9 = 7-2^2- (3-3) 



Zn, = I'd^'M exp(-Kn(M)). (3.4) 

A very useful Schwinger-Dyson identity for this pure matrix model can be easily derived 
from the path integral 

i:/''""^^E4(^'v-'"')^"■ *^-^' 

We obtain (with G^ =< TvnM'^ >„ /N^, Gf =< TtnEM'^ >m /N^ and Gf =< 

1 = G^^ + 2m2ArG^ + 45G^. (3.6) 

We first consider the case r^ = or equivalently m} = 0. In this limit the model reduces 
to the famous hermitian quartic matrix model. The partition function and the potential 
are given in this case by 



Zo = fd'^'M exp(-l/o(M)) , Vo{M) = TtnC^M'' + J^mA. 



(3.7) 



This is studied originally in the seminal work [13j. The physics of this model is character- 
ized by the celebrated one-cut to two-cut transition [H] which is expected and is known 
to persist at the level of noncommutative scalar phi-four theory. Indeed we can convince 
ourselves from Monte Carlo simulations that the one-cut to two-cut transition is lifted 
at the level of the noncommutative theory to the transition from disorder to the stipe 
phase [15]. 

Our primary reason in this article for applying the renormalization group idea to non- 
commutative phi-four is to probe this transition by means of an analytic non-perturbative 



method. We believe that the one-cut to two-cut transition can be accessed in a satisfactory 
way with the scheme put forward in this article which consists in formulating noncommu- 
tative phi-four model in the matrix basis and employing the approximate renormalization 
group recursion formula and 1/A^ expansion. All this is based on the equivalence con- 
jectured for example in [5] between noncommutative field theories at = oo and large 
N matrix models which is assumed to hold beyond the critical point. For an alternative 
viewpoint we refer to [9] who used a different regularization of noncommutative field the- 
ory based on the twisted reduced model. In any case in the disorder phase where both the 
symmetry M — > —M and translational symmetry are respected large noncommutative 
field theories are certainly equivalent to large A^ matrix models. 

For completeness we also mention here that the phase structure of noncommutative 
phi- four is far more complicated |15pi6j. Besides the one-cut to two-cut transition there 
is the usual Ising transition and also a transition from a non-uniform order to a uniform 
order which all meet at a triple point. As opposed to the usual uniform ordered phase 
where the symmetry M — > —M is spontaneously broken and translational symmetry is 
respected in the non-uniform ordered (stripe, two-cut) phase the symmetry M — > —M 
is maintained whereas translational symmetry is spontaneously broken. 

The Schwinger-Dyson equation ()3.6p for r^ = can be put into the form (with G2 = G2 
and Gl = G4) 

r2 = l + 8<7G2-4<7r4(G2)^ (3.8) 

The connected 2— point and 4— point functions and the 2— point and 4— point proper ver- 
tices for this model are given by 

C2 = G2 , Fa = {C2)-^ , G2 = ^a2(4 - a^) = 1 - 850 + 144^2 _ MhQgl + .... (3.9) 

C4 = G4 - 2(G2)^ , F4 = -C74(C2)-^ , Ga = a^(3 - a^). (3.10) 

Next we consider the case of small r^. The partition function is given in this case by 

Zm = I d^^M exp(-yo(M) - m^TrNEM^). (3.11) 

Now we will employ the method elaborated in [T7j to calculate approximately this partition 
function. Omitting the somewhat lengthy detail we obtain an effective multi trace matrix 
model (correct up to order rn^) given by 

/-. 2 

d^'M exp ( - 2(1 + 2Nm^)TrNM^ - j^Tr^M^ + '^{TrNMf)(2. 12) 

We compute immediately 

< TtnEM^ >m = -d^2 In Z„, 

= N <TrNM^>^-]^<{TrNMf>m. (3.13) 
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After a straightforward tedious calculation (see next section for a similar calculation) we 
arrive at the 2— point function (with g = g/{l + m?N)'^, h = 1 + rr?N jl + ...) 

,2 



< Mtj > 



6^JN 



1 + m'^N/2 



G2{9,h) 



m 



1 + m?N/2 



[i-l/2){l + hdh)G2{g,h) + 0{m'') 



(3.14) 



G2{g, h) = l-8gh + g'^{64:h + SO/i^) + ... 
From this result we conclude that 



G5 



1 



1 + m'^N/2 



G2ig,h) + ... 



(3.15) 



(3.16) 



g: 



2 1 + m'^N/2 



G2{9,h) + 



(3.17) 



Thus at the leading order in tti^ we have G™ = G™/2 and as a consequence the Schwinger- 
Dyson identity becomes 1 = (1 + rn?'N)G'^ + AgG^ . This also means in particular that 
at the leading order in m? we have < Trj\fj'^ >m= + 0{m?) where J = M — TtnM/N . 
The partition function can then be reduced further to 



jAT^ 



1 



d''"M exp ( - -(1 + Nm'^)TrNM^ 



N 



TtnM^ 



m^N 



TtnJ^) 



f d^'M exp ( - ^(1 + Nm^)TrNM^ - j-TrNM^). 



(3.18) 



Furthermore after a field scaling M — > M = (1 + Nm? /2)M we obtain G™ = (1 — 
Nm?)G2{g, Nm?) and G™ = (1 — 2N m?)G ^{g , N m?) . The primed correlation functions 
are computed with the partition function 

'd^V exp ( - \Tr^M'^ - J^Tt^M'' - H^Trr^j''). (3.19) 



Zn 



The Schwinger-Dyson identity becomes 1 = ^2(5, Nm'^)+'igG^{g, Nm?). In this form this 
identity holds true in an obvious way for m? = since ^2(5,0) = G2{g) and ^4(^,0) = 
G^{g)- We define the proper 2— point vertex r2(^, Nm?) by T2{g, Nm?) = G^^{g., Nm?). 
A straightforward calculation gives 



T'2{g,Nm^) = 1 + AgG'^{g,Nm')G'2~\g,Nm') 

= r2(5,0) +45iVm2a^„2(GlG2-^)|iVm2=o + - 

= r;(^,o) + iVm2(2(r; - 1) + a^™2r;u„2=o) + 

We rewrite this as 

(1 + Nm^){T'2 {g, Nm^) - 1) = (1 + iVm^) (Fs (5, 0) - 1) + Nm^ 



(3.20) 



r2-i 



(3.21) 



Let us recall that T2{g-,Nm?) — 1 is the quantum correction to the mass of the field 
M which is equal classically to exactly f. Now we undo the rescaling M — > M = 
(1 + Nm? /2)M. The quantum correction to the mass of the field M is V2{g^ Nm?) — 1 = 
(1 + Nm'^){T2{g, Nm?) — 1) while the classical mass is 1 + Nm'^ as is obvious from (|3.18|) . 
We obtain then 



T2{g,Nm'^) = T2{g,0)+Nm^ 



T2ig,0) + Nm'^ 



r2 - 1 + 9jVm2r2|Arm2 = 



r2 - 1 + ^9/^^2U=l 



(3.22) 



In the following section we will give, as a by product of our analysis of the noncommutative 
phi-four theory with oj = 0, a, direct derivation of this formula together with an explicit 
perturbative expansion for dfyfjy^2T2{g, Nm?). In fact the above formulas will be used 
to test our method. We will also see in the following that the combination Te{g,0) = 
(1 — dj\f^2T2{g,rn?))\m'2=o is precisely the proper vertex associated with the operator 
TrNEM^. 

For a systematic more sophisticated study of the Penner matrix model (j3.2p we refer 
the reader to the recent article 1191. 



4 Wilson Renormalization Group 
4.1 Recursion Formula 

In this subsection we review very briefly the original Wilson renormalization group 
recursion formula. By using the principle of truncation (more on this below) the action 
at any renormalization group step is taken to be given by 

Sn{(p) = \j d''x{d^<l)nf + f d^xPn{(t>n). (4.1) 

First we divide momenta logarithmically as 1/2 < \k\/K < 1/2 where A is the cutoff. 
Next by integrating out field modes with momenta in the highest shell corresponding to 
/ = and neglecting momentum dependence within each cell in position space (among 
other things) we arrive at an action of the same form with the replacements (pn — > (pn+i 
and Pn — > Pn+i where 

Q„+i(2'^/2a-iz) ^ ^,-ip„+i(0„+,) = -2'^ln^. (4.2) 

The variable z is related to the field (pn+i in a particular way which does not interest us 
in this article whereas 'w~^ is the volume of a single cell in position space and a„ is the 
normalization of the field. Indeed the fields (pn and (pn+i are related by 



fn 



[x) = $„(x) + 2-'^/^an^n+i{x/2). (4.3) 



The background field <^„ contains precisely the integrated momenta 1 < |fc|/A < 2. The 
function Iniz) appearing in the recursion formula (I4.ip is given by the integral 

Iniz) = dy exp(-y^- -Qn(y + z) - -Qn{-y + z)\. (4.4) 

4.2 Integration and Truncation 

The Wilson renormalization group idea consists in general in the following three main 
steps: 

• Integration and Truncation: Define a new (renormalized) action S by integrat- 
ing out all modes with momenta larger or equal than pA. where < p < 1 and A 
is the cutoff. This is a very complicated step which requires in practice the use of 
truncation and many other approximations before it can be carried out explicitly. 
Wilson truncation means that we calculate quantum corrections to only those terms 
which appear in the starting action. The principal approximation we will further use 
consists in the above Wilson recursion formula (|4.ip . We state without proof (see [1] 
for a derivation) that the use of this recursion formula is completely equivalent to 
the use in perturbation theory of the Polyakov- Wilson rules given by the following 
two approximations: 

— We replace every internal propagator 1/(A;^ + y?) by l/(c^ + //^) where c is a 
constant usually taken to be A. 

— We replace every momentum loop integral J ^dk/ {2Tr) by another constant 
vd = A ?)/) where 

. _ 2(1 - p^) 1 1 

"^ ~ D (47r)^/2 T{D/2) ' ^ ' 

We remark here that the above two rules are rather obvious and in fact natural 
in the limit of the dilatation parameter given by p — > 1 regardless of the Wilson 
recursion formula itself. 

• Rescaling: This a simple step which consists in restoring the cutoff to the original 
value A which can be achieved by simply making the change of variables p — > p = 
pp. 

• Normalization: Restoring the standard normalization of the kinetic term which is 
equal to 1/2 by appropriately rescaling the field variable. This is also a very simple 
step. 

Let us then start by decomposing the N x N matrix M into an A^ x A^ background matrix 
M and an A^ x A^ fluctuation matrix m, viz 

M = M + m. (4.6) 

The background M contains slow modes, i.e. modes with momenta less or equal than pA 
while the fluctuation m contains fast modes, i.e. modes with momenta larger than pA 
where < p < 1. By integrating out the field m we obtain the partition function 

9 



By using momentum conservation, the symmetry m — > —m and also by neglecting non- 
planar diagrams (in the sense of the 1/N expansion) we find that the non-perturbative 
correction AS'(M) is given explicitly (up to the fourth power in the field M) by the 
cumulant expansion 

A5(M) = 4^ f d^x<TrM^m'^{x)>co 

- 8-^ / d^x / d^y < TrMm^{x)TrMm^{y) >co 

- 8^ I d^x I d^y < M^m^{x)TrM'^m^{y) >co 

+ 32^ / d^x / d^y / d^z < TrMm^{x)TrMm^{y)TrM^m^{z) >co 

- — ^ / d'^x / d'^y / d'^z / d^w < TrMm^{x)TrMrr?{y)TrMm^{z)TrMm^{w) >co • 

(4.8) 

The notation "co" stands for the connected component. The first and second terms 
yield correction to the mass parameter whereas the last three terms yield correction to 
the quartic coupling constant. The wave function renormalization is obtained from the 
expansion around p^ = of the second term. This is the most difficult contribution to 
calculate as we will see shortly. 

Mass and Harmonic Oscillator Renormalizations: Quantum corrections to 
the mass parameter ^^ and to the harmonic oscillator coupling constant r^ are obtained 
from the first term of (j4.8p and also from the second term of ()4.8p evaluated at p^ = 0, 
viz 

A5„,ass+h.o = -^ d^x< TrM^m^x) >co -^ / d^x / d^y < TrMm^{x)TrMm^{y) >,o |p2=o- 

(4.9) 

The corresponding Feynman diagrams are shown on figures ([6]) and ([7|. We included 
diagrams up to three loops. The series can be continued if one wishes to go further. 
However all the subleading (in A) Feynman diagrams which dominate the large N limit 
can be recovered in our scheme as we will now show. By using first the approximations 
(j2.12p and (j2.13p and then the Wilson recursion formula as described above the first set 
of Feynman diagrams shown on figure ([6]) reduces to 



4^ f d^x<TrM^m^{x)>co = f d^xM^{x)i 



4vDuAnai{c) - 32{vdu) A„yj(c) A„o„q (c 



(4.10) 
10 



By expanding the second set of Feynnian diagrams shown on figure ([7]) around p^ = and 
retaining the th order term we get extra corrections to the mass and harmonic oscillator. 
By using again the approximations (12.121) and (I2.13P and the Wilson recursion formula 
we get explicitly 



-8^ / d^^ / d^y < TrMm^{x)TrMm^{y) >co |p2=o = / d^xM{x)ijM{:, 



(4.11) 

The terms multiplying M'^{x)ii and M{x)ijM{x)ji still depend on the indices i and j 
which is the source of the harmonic oscillator renormalization. Recall that the harmonic 
oscillator term is of the form J d^xM'^{x)ii{i — 1/2). The mass correction corresponds to 
setting these indices equal to some fixed value n whereas the harmonic oscillator correction 
corresponds to expanding these terms linearly around n. Again in the spirit of Wilson 
contraction we treat the index i (or j) as a continuous variable and expand the propagator 
around i = n as follows 



Ainoic) = A„„„(c)-(i-n)r2AL„(c) + (i-n)VALo(c) + ... (4.12) 

This actually makes sense since we are assuming that is sufficiently large and thus r^ is 
sufficiently small. 

The full mass and harmonic oscillator renormalizations are contained in the sum of 
the two terms (j4.10p and (j4.1ip . By using the expansion (j4.12p we obtain after some 
calculation the result 



^'-^ma.ss-l-h.o — I U X 



^(c2 + fi^ + r^N)r2{g, h)TrNM\x) + K^d^r2{g, h)TrN{E -n + ]^)M\x) 



n r 



4 



2 c2 + ^2 



diT2{g,h)TrNEM\x) 



+ ... (4.13) 



This formula contains the leading as well as the subleading corrections in 1/9 to the mass 
and to the harmonic oscillator terms. The 2— point proper vertex r2(g, h) is defined by 

T2{g, h) = 8gh - SOg'^h^ + Ig^ilbGh^ + 576/i^) + .... (4.14) 

The new coupling constants g and h = 1/k, are defined by 



,2 /„\ u — ^nno{c) 
AnQno\C) 



g = v^uAl^^^ic),h = ^p^^. (4.15) 



At (9 = oo we have /i = 1 and thus T2{g, 1) = 8g - 80g'^ + 1664/ + The 2-point 

proper vertex T2{g, 1) is expected to be related to the 2— point proper vertex T2{g) of the 
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quartic matrix model, i.e. the expansion 8g — 80g^ + 1664g^ + .... should be recognized 
from the pure quartic matrix model. Indeed we have 

T2ig)-l = T2ig,l). (4.16) 

The 2— point proper vertex r2(5') is in fact a function known non perturbatively given by 
(with a = (VI + 485 - l)/245) 

r2(<7) = ^2u_ a2) = 1 + % - 8O5' + I66453 - .... (4.17) 

The full mass and harmonic oscillator action is therefore given by (using also h = 1 + 

5mass+h.o + A5^ass+h.o = ^ (^' + 0^)^2(9) " c^ + r^N{r2{g) - Te{g)) + 0{r^)\ j d^'xTrNMH^ 

+ r^heig) - 0^^^-re(5,/i)U=i + 0{r^)\ j d^xTrNEM\x). 



(4.18) 



We have defined 



Te{g,h) = i + ^a,r2(ff,/i) , r,(5) = re(5,i). (4.i9) 



In contrast to T2{g) the function Tg^{g) is only known perturbatively. Note also that we 
did not need to fix the value of n in deriving the result ()4.18p . In the following we will see 
more clearly the virtue of the choice n = 1/2 which was assumed in a previous section on 
pure matrix models. For c^ = the mass renormalization agrees with the result obtained 
from the Penner matrix model. 

Wave Function Renormalization: The wave function renormahzation is also ob- 
tained from the 2nd term of ()4.8p and as a consequence the relevant Feynman diagrams 
are still given by those shown on figure ([7]). More precisely we need, as before, to expand 
these diagrams around p^ = but retain now the linear term in p^ which is very difficult 
to do explicitly. We will present here a new estimation of the coefficient of p^ which is 
consistent with the recursion formula. 

The three graphs shown explicitly on figure (jT]) are of the form 

I ^ J ^^^-o{P+Pi+P2)f{pi,P2). (4.20) 

We need the expansion around p^ = of the propagator Aj„f,(p + pi +P2) = l/((pi +P2 + 
p)"^ + c'^) where a'^ = fi'^ + r^ {i + uq — 1) . By using rotational invariance we can replace 
^inoip + Pi +P2) in this integral as follows (with e = 4 — D and eo = e/D) 



^inoiP + Pi + P2) = A,„„(Pi +P2) + [eo{pi+P2f - <t^]aL„(pi +P2)/ + ■•■• (4.21) 
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The first term corresponds to the p^ = contribution relevant for mass and harmonic 
oscillator renormalizations whereas the second term will provide wave function renornial- 
ization. The contribution of this last term to the integral (I4.20p is given by 



P 



I ^ I ^Hpi+P2f -a']Al^{p,+P2)f{pi,P2). (4.22) 



In the recursion formula we replace the internal momentum appearing in each propagator 
by an average value c which we have chosen to be c = A. We anticipate the application 
of the recursion formula and replace the factor eo(pi +^2)^ — <7^ with eocP — a"^ . We get 
then the estimation 



P^ [eoc 



a 



d pi f d^p' 



{2Tr)D J (27r) 



^fnniPl+P2)fiPl,P2)- 



\£) "ino 



(4.23) 



In higher order diagrams there could be more than one propagator depending on p and 
also the propagators may depend on more than just two internal momenta which will 
complicate the situation even further. More precisely if we generalize ()4.2ip to higher 
order terms than the extra factor eo(pi +P2)'^ — c^ may be different. This important issue 
will be discussed further below. 

The graphs shown on figure ([7]) then yields the result (using also the approximations 
(j2.12p and (|2.13p and the recursion formula) 



i^ I d^x I d^y < TrMm^{x)TrMm?{y) >c 



(2^ 



p'^Mij{p)Mji{-p) 



fi^ + r^{i + no-l) 



eoc 



VDu)'^Aljc)Anoj{c)Anonoic) " 192ivDU 



X A^ ^(c)A2 (c)A2 .(c) - U{vDufAljc)AiJc)^^ 



^2 ^ 



c +.... 



(4.24) 



We split the factor p^[//^ + r'^{i + uq - 1) - eoc^] as p^[/i^ + r'^{no - 1/2) - eoc^] plus 
p^r'^ii — 1/2). The first bit depends on the indices i and j only through the propagators. 
Now we only need the contribution obtained by setting these indices equal to the fixed 
value n = 1/2. The linear term around n = 1/2 will lead to a correction of the form 
Tr]\fE(dfj_M)'^ which is not present in the original action. For n 7^ 1/2 the correction 
will be of the form Tr^^E + constant) (9^M)2. By the same reasoning we discard the 
correction Mij{p)Mji{—p)ijP'r'^{i — 1/2)) since it corresponds to a correction of the form 
Tr]\[E(dfj_M)'^. This is the principle of Wilson truncation in action again. We get then 
the result 

"^""^ p'Mij{p)M,i{-p) 



8^ / d^x / d^y < TrMm^{x)TrMm?{y) >^ 



(2vr) 



D^ 



X [/x^ + r^(no - -) - eoc']\ 8(t;Dn)^A^^„(c)A„„„„(c) - 256(t;Dn)^A^„„(c)A^^„^ic) 



+ 



(4.25) 
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We rewrite this as 

- 8^ / d""- 1 d^y < TrMm\x)TrMm%y) >eo= ^l^f,^i)~^"^2 [^^(5, 1) 

dhZ2{9,h)\h=i+0{r^)] j d^'xTrNid^Mf. (4.26) 



5 N 



2(c2 + ^2) 
We have defined 

Z2(ff, h) = Sg'^h^ - 256/ /i^ + .... (4.27) 

We need now to find a combination of Green's functions and proper vertices of the pure ma- 
trix model with an expansion given exactly by 8g'^ — 256g'^ + .... Prom the Schwinger-Dyson 
equation (j3.8p we propose that the function 2gT4{g)G2{g) is the correct guess. Notice the 
resemblance of the graphs corresponding to 85^2 and —4:gr4G2 to the graphs associ- 
ated with the terms 4:{u/N) f d^x < TrM^m^{x) >co and -8{u'^/N'^) J d^x J d^y < 
Tr Mm^ {x)TrMm^ (y) >co respectively. Indeed we compute 

^2(5) = ^2(5,1) = 2gU{g)GUg) 

2ga'^{l-a'^){5-2a'^) 



3 4-a2 

Q^2 oc;c„3 



256ff^ + ... (4.28) 



The full kinetic action is given by 



•Skinetic + AS'kinetic = S ^T H 2 i "2 -^^(g) + „. 2 1 2'|2 ('^^(^ + ^o)Z2{g) + {iJ''^ - CqC^) 

X dhZ2{g,h)\h=i)+0{r^)\ I d^'xTrNidf.Mf. (4.29) 

As mentioned previously in higher order diagrams there could be more than one propagator 
depending on p and also these propagators may depend on more than just two internal 
momenta and as a consequence the extra factor eo(pi+P2)^— o"^ appearing in (I4.23p may be 
different. This coefficient leads to the crucial factor (/i^ — eoc^)/(c^ + fJ-'^) which multiplies 
the non perturbative function Z2{g). The central conjecture of this section is that within 
the scheme of the recursion formula outlined above the coefficient (//^ — eoc^)/(c^ + ^J^) is 
in fact the same coefficient appearing in higher order corrections of the function ^2(5, h). 
In order to introduce a non trivial wave function renormalization we can alternatively 
employ the simpler approximation motivated by dimensional consideration used in the 
context of matrix models in [lOpilj . In this scheme the first derivative of the propagator 
with respect to the external momentum p^ will be approximated by the multiplication 
with the given propagator as follows 

P^T^(---)lp2=o = -P^T — , x2 , 2 (---)lp^=o- (^-30) 

ap^ '^ ypi +P2+ p) + "^ 

In this case the integral (j4.20p is approximation with 

"^ y W^ J (^^»o(Pi+:P2)/(pi,P2). (4.31) 
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The difference between (|4.23p and (j4.3ip consists in the difference between the factors 



Finally we also mention here Golner extension of the recursion formula derived in |12] 
which also allows a non trivial wave function renormalization. As it turns out this is 
a more consistent framework which requires however a non trivial modification of the 
starting action and thus we do not discuss it any further. 

Quartic Interaction Renormalization: This is obtained from the last three terms 
of equation ()4.8p and without any further calculation we simply state the final result 
obtained using the approximation (j2.12p and (j2.13p and the recursion formula. This is 
given by 



ASinteraction = -[^^^ - 1] j d"" xTtnM^ (4.33) 



The function T4{g, h) is now defined by the expansion 
'^4(9, h) 



45 
The full interaction term is therefore 



1 - 8gh^ + 5^(160/1* + 6Ah-') + .... (4.34) 



'interaction ~r ^"Jintcraction — ,t . I "• •^-' ^-'" (X 



N 4g J 

|j^{r4(g,l) + ^(/^^,^ 9.r4(g,/.)U=i + 0(r^)}|d^xTrM^C 



X). 

(4.35) 



As expected the 4— point proper vertex T4{g, 1) is precisely the 4— point proper vertex of 
the quartic matrix model given by the expansion 

. , , , 9(l-a2)(5-2a2) ■? ^ . x 

U{g, 1 = U{g) = ^ ^,/^ ,.^ ' =Ag- 32^^ + 8965' + ... 4.36 

4.3 Renormalization Group Equations 

So far we have only gone through the first step of the Wilson renormalization group 
approach, namely integration and truncation. By putting the results (I4.18J) . (J4.29P and 
(j4.35p we obtain the effective action including all leading and subleading corrections in 1/6 
at the leading order in 1/A^. By performing the second step of the Wilson renormalization 
group approach, i.e. by scaling momenta as p — > p/ p so that the cutoff returns to its 
original value A we obtain the effective action 



15 



S + AS 



2 
rV|re(9) + 



r^Af 



(c2 + /i2) 



2\2 



AZ(5,M') + 0(r4) 



(27r)^ 



p^TrNM*{pp)M{pp) 



T{g,fi^)+r^NAT{g)+0{r'^ 



(2^ 



TrNM*{pp)M{pp) 



r'^{N +1) 



up 



^A{g) 



N 4g 



2-^re(5)+0(r4) 
-^^AU{g) + 0{r' 



A jD 



d^p 



(2vr) 



^TrNEM*{pp)M{pp) 
1 /^rf^P4 







(2^)^-7o (2vr) 



\D 



TrNM{ppi 



The functions Z{g,fj?), r((7, ^2), r2((7) and r4(g') are known non perturbatively whereas 
we were able to determine the functions AZ{g,p'^), AT{g), Ar4{g), AFeig) and r^ig) 
only perturbatively. These functions are summarized in table ([1]). 

The third and final step of the Wilson renormalization group approach consists in 
rescaling the field in such a way that the kinetic term is brought to its canonical form as 
follows 



M (p) 



2+D 
p 2 



^Z{g,p^) + 



r^N 



(c2 + H 



2^2 



AZig,p^)M{pp). 



(4.38) 



We get the action to the form 



S + AS 



1 



2 70 (2^)^ 



"^ ^ p^TrNM'*{p)M\p) 



1 

+ ^ 



P 



2Zig,p^) 
^ d^p 



T{g,p^) + r^NAr{g) 



+ 



(2vr)^ 

p-2^2 



Z{9,P^ 

A jD 



TrNM*{p)M (p) 
r^{N+l) 



Teig) + 



C^ + p^ 



Are(5) 



r^N Tig,p^)AZig,p^) 
(c2+^2)2 zig,p^) 

r^N T,{g)AZ{g,p^ 



+ 0{r^ 



d^p 

W)^ 



TrNEM*{p)M (p) 



+ 



u 1 



P 



NAgZ^{g,p^ 



r^ig) + 



(? + /i2 



Ar4(ff) 



(c2+/.2)2 Z(<7,/i2) 

2r2Ar r4((7)AZ(5,^2 



+ O(r^) 



A aD 



d'^p^ 
(2vr)^-yo (^^ 



(c2+/.2)2 ^(^^^2 



Tr^vM (pi)...M (P4)(2^)^<5^(pi + ...^4). 



+ 0(r^ 
(4.39) 



Let us recall that r^ is related to the harmonic oscillator coupling constant Q by 2^2^^ = 
47r(ri2 + 1) and that we have started with an action at the self-dual point, i.e. 0,^ = 1 
or equivalently uj = 0. Thus the renormalization of r"^ can be understood either as a 
renormalization of ^2 or as a renormalization of the noncommutativity parameter. Both 
descriptions must be equivalent but we prefer to think that we have always Q^ = 1 and 
attribute all the renormalization to 6. 



■M{ppi) 
(4.37) 
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In position space the above action reads as follows 



S + AS = )- f d'^xTrNid^M'f + ^ / d^xTr^M'^ + r'^ / d^xTrNEU"^ 



+ ^ / d^xM'4 



(4.40) 



The renormalized mass ^ ^, the renormalized quartic coupling constant u and the renor- 
malized inverse noncommutativity r ^ are given by 



At 



'2 _ P'^ 



Z{g,lJ? 



T{g,^?)+r^NAT{g) 



r^N Tig,fi^)AZig,fi^ 
(c2+^2)2 z{g,^^^) 



+ 0(r'' 



■ (4.41) 



'2 



p ^r^ 

Z{g,iJ?, 



^e{g) + ^^^^^^^Are(5) 



r^AT re(ff)AZ(<7,^2 



(c2 +^2)2 



Z{g,fi^ 



+ 0{r^ 



. (4.42) 



u = u 



Z^{g,p?)Ag 



T^{g) + ^— ^Ar4(<7) 



2r2iV r4(5)AZ(<7,^2 



(c2 + ^2)2 



Z[g,fJ-^ 



+ 0{r^ 



(4.43) 



The process which led from the bare coupling constants /U^, r'^ and n to the renormalized 
coupling constants fj,"^, r'^ and u can be repeated an arbitrary number of times. The bare 
coupling constants will be denoted by /^q, Tq and uq whereas the the renormalized coupling 
constants at the first step of the renormalization group procedure will be denoted by fif, 
rf and ui . At a generic step / + 1 of the renormalization group process the renormalized 
coupling constants ^f^,^' '"f+i 



and n^+i are related to their previous values iJ.f, rf and ui 



by the renormalization group equations 



P'f+i 



Zigup-i] 



Tigi,fif)+rfNATigi 



rfN 



r{gi,fif)AZ{gi,fil 



(c2 + nf 



2n2 



z{gi,pf) 



(4.44) 



r2 






re(5z) + 



rf(iV+l) 



2 2 -^re(<7z) 

C2 + pf 



rfN 



T,{gi)AZ{gi,pl 



{c^ + pfr 



ZigufJ-l 



(4.45) 



Ul+l = ui 



Z^gi,pf)^g 



Uigi) + 



rfN 



^T^{gi) 



2rfN U{gi)AZ{gi,pf] 



(c2 + ^2)2 



Z{gi,pf) 



. (4.46) 



The effective coupling constant gi will of course be given in terms of pf, rf and ui by 
the same formula that related ^o to p^, Tq and uq. We are therefore interested in renor- 
malization group flow in a 3— dimensional parameter space generated by the mass p'^, 
the quartic coupling constant u and the harmonic oscillator coupling constant (inverse 
noncommutativity) r^. 



17 



5 Fixed Points and Critical Exponents 

By definition a renormalization group fixed point is a point in the space parameter 
which is invariant under the renormalization group flow. If we denote the fixed point by 
fil, r1 and u^, then we must have 



A^* 



P 



Z{g*,iJ.l) _ 



ng,,^ll) + rlN^V{g,) 



rlN T{g,,^ll)^Z{g,,^ll 



(c^ + ^i^. 



2^2 



^(5*,M*) 



(5.1) 



Z{g*,nl) 



r (n)^ !:M±}lAr (n ^ - "^'^ re(g.)AZ(ff,,/i^) 

^^'*'^ c' + pI ^^'*^ {c' + pI? Z{g.,^^l) 



(5.2) 



u^ = u 



*Z'^{g^,fil)Ag^ 



r4(5'*) + 



c2 +^2 



Ar 



A{g* 



2rlN T^{g,)^Z{g,,^ll 
(c2+^2)2 z{g,,i,l) 



. (5.3) 



The second equation is new by comparison with the commutative theory. The definition 
of g^ in terms of /i^, r^ and n^, is obvious. There are possibly several soultions (fixed 
points) of interest to these renormalization group equations. We will mainly concentrate 
on the matrix model fixed point corssponding to infinite noncommutativity which is the 
most obvious solution to equation (|5.2p given by 



rl = 0. 



(5.4) 



The remaining two equations reduce then to 



P* 



Z{g*,pl 



:^{g*,pl) 



(5.5) 



u* = u* 



Z'^{g*,pl 



-^i{g*)- 



(5.6) 



Thus this fixed point is fully determined by functions which are known non-perturbatively. 
An obvious solution to (15. 6|) is li* = which corresponds to the usual Gaussian fixed point. 
By discarding this solution equation (j5.6|) becomes 



1 



P 



^i{g*) 



Z'^{g*,pl) ^g* 



We define the new variable 



(5.7) 



a 



P 



(? + ^2 



The solutions of equations (15. 5p and (|5.7p are given respectively by 

fi{p,g*) + 2eoZ2{g*) 



a* 



4(l + eo)^2(5*) 



(5.8) 



(5.9) 
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4(l + eo)Z2(<7* 
The functions /i and /2 are defined by 



«. = hiP^9*) + ^loM9*) _ (5_^o) 



/i = -1 + ^ - a/(1 - ^ - 2eoZ2(r7*))2 + ^(1 + eo)Z2ig.)(T2ig.) - 1). 
p^ V p p 



(5.11) 



/2 = -2 - 2^-^^^ + 2Z2(5*)6o. (5.12) 

Clearly we must have fi = f2- In other words g* is the solution of the consistency equation 

f = fi-f2 = 0. (5.13) 

The physical region of g is [0,oo[ while the full domain of definition is [— 1/48, oo[. Fur- 
thermore the functions r2, r4 and Z2 depend on g only through a = a{g) defined by 
a^ = 2/(yT+^8^ + 1). Graphically we observe that the two functions /i and /2 intersect 
for all dimensions d = 2, 3, 4. We will employ the Newton-Raphson algorithm to determine 
the roots of the consistency condition (j5.13p . The roots are then given by the iterative 
equation (with an initial step which in principle can be anything) 

After a sufficient number of iterations n we will obtain the root Xn = g*- We determine a^, , 
the critical mass fil = At^/A^ and the critical coupling u^, = u^K^~^ from the equations 

/2 + 2eoZ2(ff„) 

pI = r^- (5.16) 

1 — a* 

u. = Mi±^. (5.17) 

VD 

The results of this calculation are shown on table ([2]). For comparison we report the critical 
values in the approximation (j4.3ip in table ([3|). There is of course in each dimension the 
extra Gaussian fixed point as we have discussed. There is only the Gaussian fixed point in 
D = A{d = 6) in the approximation (j4.3ip which can not be true since the matrix model 
fixed point is expected to exist in all dimensions. In our approximation (|4.23p we have 
checked that there is always a non-trivial fixed point for any value of p in the interval 
0<p< 1. 

In the remainder we will compute the mass critical exponent v and the anomalous 
dimension r/ within this scheme. 
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Theoretical Digression: The computation of the mass critical exponent v requires 
hnearization of the renormahzation group equations (14.441) , (I4.45P and (I4.46P . Let us first 
rewrite these renormahzation group equations in the compact form 

G("+i)=7W(G("\p). (5.18) 

The vector of couphng constants G is defined by G = (^1,^2,^3) where G\ = jj.'^, 
G2 = u and G3 = r^. The hnearized renormahzation group equations are of the form 
(with 6G = G-G^) 

6G^''+^'> =M{G,,p)6G'^''l (5.19) 

The matrices M and A4 are different for different fixed points. The matrix M is related 
to M as 

M^,{G,,p) = ^A^,(G("),p)|c(„)=G,- (5.20) 

Let fa (p) be the left eigenvectors of M with eigenvalues Xaip), viz 

vl(p)M{G,,p) = K{p)M{G,,p). (5.21) 

The so-called scaling fields Ua are defined by the projection of the vector of coupling 
constants 5G onto the left eigenvectors v'^{p), viz 

Mp) = vl{p)5G. (5.22) 

These vectors do not mix under the renormahzation group transformations since they 
satisfy 

u^:+'\p) = K{p)u^:\p). (5.23) 

Thus the scaling field Ua will increase under renormahzation group transformations (rel- 
evant) if Aq > 1 while it will decrease if Aq, < 1 (irrelevant). 

It turns out that every real asymmetric matrix M can be factored (Schur factorization) 
as M = U^^TU where T is upper triangular and U is unitary. This substitutes for 
diagonalization which is not always possible for real asymmetric matrices. The eigenvalues 
in this case are real or come in complex conjugate pairs. As usual the eigenvalues are the 
zeros of the characteristic polynomial. 

A fundamental property of the renormahzation group is the semigroup property which 
for two different dilatation parameters pi and p2 reads 

M{G„pi)M{G„p2) = M{G„p2)M{G,,pi) = M{G,,pip2). (5.24) 

It is not difficult to convince ourselves that this leads to the requirement that the renor- 
mahzation group eigenvalues Xa{p) scale as 

A«(/9) = A,(l)p-^". (5.25) 
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The so-called scaling indices j/q, (also called critical exponents) are the renormalization 
group eigenvalues of the scaling fields Ua{p)- Indeed under a renormalization group 
transformation we have Ua{p) — > u^{p) = \a{p)ua{p) and thus for an infinitesimal 
renormalization group transformation {p — > 1) we must have u^{p) — \a{l)ua{p) = 
—ya^"^ pXaO-)ua{p) + 0((lnp)^). In other words 



P-r^a{p) = -VaUaip). (5.26) 

dp 

In the limit p — > we have Ua{p) — > ±00 when ya > (relevant) and Ua{p) — > ±00 
when Ha < (irrelevant). 

The Mass Critical Exponent z/: The matrix M in our case is of the form 

Mil M12 Mi3 

M21 M22 M23 I . (5.27) 

M33 

The corresponding characteristic polynomial gives the two equations 

The eigenvalue in the direction G3 = r^ is therefore given by 

This eigenvalue (or more precisely In A3) is plotted on figure ([1]) as a function of In p. It 
is immediately obvious that A3 > 1 and hence r^ is a relevant coupling constant like the 
mass. However the function Te{g) used in the above formula is only known perturbatively 
and hence this conclusion should be taken with care. 

The two remaining eigenvalues are determined from the linearized renormalization 
group equations in the 2— dimensional space generated by Gi = p? and G2 = u. As it turns 
out this problem depends only on functions which are fully known non perturbatively. The 
eigenvalues Ai and A2 can be determined from the trace and determinant which are given 
by 

Ai + A2 = Mil + M22 = Tr2M , A1A2 = M11M22 - M12M21 = det2M. (5.30) 

In other words 

^^ ^ Tr.M ± V(T,.M)^ - 4det.M _ ^^ ^ ^^^^ _ ^^^ ^^^^^ 

These are real as long as (Tr2M)^ — 4det2M > 0. Indeed a real asymmetric matrix can 
have in general complex eigenvalues. For a discussion of the relevance and interpretation 
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of complex renormalization group eigenvalues see chapter 3 of [20] and references therein. 
This point is not discussed in this article. 

The formula (I5.25P was used as a crucial test for our numerical calculations. In partic- 
ular we have determined by means of this formula the range of the dilatation parameter p 
over which the logarithm of the eigenvalues scale linearly with In p. It is natural to expect 
this behavior to hold only if the renormalization group steps are sufficiently small so not 
to alter drastically the infrared physics of the problem. 

Some results of this test are shown on figure ([2]). In the approximation (j4.23p the 
eigenvalues Ai and A2 become complex and conjugate to each other at around Inp ~ — 1 
and thus the regime Inp > — 1 is of no interest to us here. In the regime Inp < —1 both 
Ai and A2 are real. The eigenvalue ln|Ai| as a function of Inp is linear in the regime 
Inp < —1 whereas the eigenvalue In IA2I is linear in this regime only for D = 2. For D = 3 
and 4 In I A2 1 is linear only for In p < < — 1 which is not very satisfactory. Fortunately it is 
the eigenvalue Ai which is associated with a relevant scaling field and thus used to define 
the mass critical exponent as we will discuss shortly. 

Note that the eigenvalues Ai and A2 are both negative in the regime Inp < —1 which 
means that the scaling indices yi and 1/2 are mildly complex. Furthermore we note that 
over the regime Inp < —1 the eigenvalue A3 scales also correctly with p (figure ([2])). 

For D = 2 the renormalization group steps can be thought of as small and the behavior 
(j5.25p holds in the regime Inp < — 1. As it happens this is the most important case 
corresponding to d = 4. We find explicitly the following fits 

In |Ai| = -1.296 Inp + 0.412 , In IA2I = 0.435 Inp + 2.438. (5.32) 

This result is shown on figure ([3]) together with the results for In |Ai| for D = 3 and D = 4 
and their numerical fits. We also obtain positive slope for ln|A2| for D = 3 and D = 4 
although the range is much more smaller. We conclude immediately that in all considered 
dimensions the scaling field ui corresponding to the mass is relevant while the scaling 
field U2 corresponding to the quartic coupling constant is irrelevant. This is the usual 
conclusion in d = 4. The critical exponents in D = 2 are given respectively by 

yi = 1.296 , y2 = -0.435. (5.33) 

The mass critical exponent is given by the inverse of the critical exponent yi associated 
with the relevant direction, viz u = l/yi- The corresponding results for yi and u are 
included in table (JH). It is very amusing to note that the results shown in table ^ are 
very close to the average value of 2/d and 2/D, viz 

. = i + i. (5.34) 

In the approximation ()4.23p we obtain instead i^ = 2/D in the limit p — > [llj . 

An alternative definition of the critical exponents which are associated with the eigen- 
values Xa is given by the equation 

Inp 

^a = - , I, I - (5.35) 

m I Aq, I 
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This is a more naive definition (to our mind) since it does not extract properly the de- 
pendence of ln|Aa| on Inp although it is the one that has been used in [TOllIT]. The 
corresponding results are shown on figure ([4]). 

For completeness we have also checked what happens in the approximation (|4.23|) . In 
this case the eigenvalues Ai and A2 are always real with Ai > A2. The eigenvalue In |Ai| as 
a function of In p is linear throughout the full regime of In p whereas In | A2 1 is linear in the 
two regimes Inp << —1 and Inp > —1 separately. We skip writing down here explicitly 
the various results for simplicity. 

The Anomalous Dimension rj: Let us go back to the wave function renormalization 
contained in the equation 



M {p) = p—Jz{g, p2) + AZ{g, fi^)M{pp). (5.36) 



We write this as 



- 2 + D-rj - I 

M{pp) = p 2— M {p). (5.37) 



The coefficient i] is called the anomalous dimension. It is given by 

ln[Z(5,p2) + ^i^AZ(5,p2)] 

r] = \ ■ 

mp 

At the fixed point we obtain 

\n[Z{g,,pl)\ 

V = ] 

mp 

ln(p-T4(g.)/4ff,) 

21np 

e ln(r4(g.)/4g,) 

2 2 In p 



(5.38) 



(5.39) 



The anomalous dimension r/ as a function of Inp in the regime Inp < —1 in the approx- 
imation ()4.23p is shown on figure ([5]). We observe that rj approaches a constant value as 
Inp — > —1. If we allow ourselves to include also the regime in which the renormalization 
group eigenvalues Ai and A2 become complex conjugate to each other then we will see 
that 7] will approach e/2 as p — > 0. This has already been observed in the approximation 

(lOTTi in nnmi]. 

6 Conclusion and Outlook 

In this article we have presented a study of phi-four theory on noncommutative spaces 
with only two noncommuting directions at the self-dual point using a combination of 
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Figure 1: The eigenvalues A3 in the approximations fl4.23p and (14.3 ip respectively. 
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Figure 2: The eigenvalues Ai and A2 in the approximations fl4.23p and (14.311) respectively. 
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approKimation (4.23) 
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Figure 3: The eigenvalues Ai and A2 in the approximations (14.23^ in the regime Inp < — 1. 
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Figure 4: The critical exponents vx and v^ in the approximation (14.231) in the regime Inp < — 1. 
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Figure 5: The critical exponent 77 in the approximation fl4.23p in the regime Inp < — 1. 
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function 


non-pert urbative 


pert urbative 


g 


Vou{c^ + li'' + r^N)-^ 




h = 1/k 


(^ + /i^ + r'^N){c^ + ^l^ + r^N/2)-^ 


l + r'^NI{2(c^ + li^)) + ... 


a' 


(Vl + 48.9-l)/249 


1- 12.9 + 288.9^ + ... 


z 


(/i' - eoC^)/{c' + /i^) 




Z2(g) = Z.2[g,l) 


((I_a2)2(5_2a2))/(18a2(4-a2)) 


8.9^ - 256.9^ + ... 


Z{g,fJ^') 


1 + 22^2(5) 




r^Cg) = r2(g, 1) + 1 


3/(a^(4-a^)) 


1 + 8.9-80.9^ + 1664.9^- ... 


r(5,M') 


(m' + c^)r2(.9) - c' 




r4(.g) = r4(.9, i) 


i9il-a'){5-2a'))/{aHA-ay) 


4.9-32.9^ + 896.9^ + ... 


Z2{9,h) 


7 


9,g'^H^ - 256g^h-^ + ... 


Ug^h) 


? 


8gh - SOg^h^ + 2.9-'{256/j.^ + 576/r'') + ... 


T4{g,h) 


7 


4.9(1 - Sgh^ + .92(160/^'' + 64/j3) + ...) 


re{g,h) 


i + 4r2(.9,/j.)/2 




Ug) 


re(9, 1) 




AZ(.g,/.2) 


c2(l + e„p')Z2{g) + i^i^ - eop'^)d,Mg, h)\,,=. 




Ar(.g) 


T2{g)-T,Ag) 




Ar4(.9) 


drXi{g,h)\h=il2 




Ar,(.9) 


-dj,{g,h)l^,l2 





Table 1: The different functions appearing in the effective action f l4.39p . 



D{d) 


9* 


M? 


M» 


2(4) 


1.845 


-0.693 


2.92 


3(5) 


2.235 


-0.732 


10.853 


4(6) 


2.518 


-0.757 


50.083 



Table 2: The critical values for p = 0.5 in the approximation f l4.23p . 



the Wilson approximate renormalization group recursion formula and the solution to the 
corresponding zero dimensional matrix model at large N. The most important results of 
this work are: 

• The Penner matrix model Tr{M'^/2 + rin?EM'^ + gM'^/N) can be systematically 
solved by the multi trace approach of [17] as illustrated by the computation of the 
2— point proper vertex in section 3. This might be an alternative approach to the 
one pursued in [19]. We hope to be able to return to this point in the future. 

• We have conjectured in this article that wave function renormalization within the 
scheme of the recursion formula is fully encoded in the approximation (j4.23p . 

• We have found a fixed point solution of the renormalization group equations (|4.44p . 
(|4.45p and (j4.46p with r^ = 0, /^^ < and n* > in ah dimensions D = 2,3,4 
corresponding to the dimensions d = 4, 5, 6 respectively. We do not exclude here the 
possibility of the existence of other fixed point solutions even within this scheme. 

• The renormalization group eigenvalue A3 was determine to be larger than 1 and 
hence r^ is a relevant coupling similar to the mass. This statement was however 
derived from the available perturbative knowledge of the proper vertex Fg and thus 
should be taken with care. 

• We have found that the renormalization group eigenvalues Ai (corresponding to the 
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D{d) 


g* 


f^l 


Ut 


2(4) 


2.282 


-0.851 


0.854 


3(5) 


0.409 


-0.643 


3.527 


4(6) 












Table 3: The critical values for p = 0.5 in the approximation (I4.3ip . 



Did) 


yi 


v 


2(4) 


1.296 


0.772 


3(5) 


1.851 


0.540 


4(6) 


2.345 


0.426 



Table 4: The critical exponents yi and u in the approximation ( I4.23p . 



mass, relevant) and A2 (corresponding to the quartic coupling, irrelevant) become 
complex conjugate of each other at the value of the dilatation parameter Inp ~ —1. 
After extracting the dependence on p we obtain the mass critical exponent 



1 1 

z/ ~ — I . 

d D 



(6.1) 



In D = 2 (the most important case for us) it is found that both the eigenvalues In Ai 
and lnA2 scales linearly with Inp over the entire real range Inp < —1. We obtain 
the indices 



yi = 1.296 (u = 0.772) , ya 



-0.435. 



(6.2) 



The anomalous dimension was found to scale as 

4- D 



r] 



p^O 



(6.3) 



We conclude this article by indicating that the most natural extension of this work should 
be to repeat the same analysis of the matrix model fixed point by replacing the recursion 
formula with the exact functional renormalization group method. The function renormal- 
ization group (as opposed to other forms of the renormalization group) is the most direct 





Figure 6: Feynman digrams contributing to the renormalization of the mass parameter and the 
harmonic oscillator coupling constant. 
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Figure 7: Feyninan digrams contributing to the wave function renormalization (linear term in 
p^) and also to the renormalization of the mass parameter and the harmonic oscillator coupling 
constant (the p^ = term). 

implementation of the Wilson idea which goes along the lines of the recursion formula but 
is exact although explicit calculation will undoubtedly involve various truncations which 
by their nature are also approximate. 
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